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Because the coefficient of thermal expansion of fused silica^ is so small (~ 0.4 x 
10~^°K~^ at room temperature), it is evident that the variation in refractive index^ 
{dn/dT ^ 1.1 X 10~^°K~^) must result from a change in the internal structure of the 
material with temperature. Note, also, that dilution of dispersive material by thermal 
expansion would tend to produce a negative dn/dT, opposite to what is observed. Here 
we develop the most elementary possible model. It is qualitatively, but not quantitatively, 
consistent with the data. 

We write the complex dielectric constant e{uj) as the sum of real and imaginary parts 

eM = e'H + ze"M, (1) 

with t'iuj) and t" {uj) each real. Then we use one of the Kramers-Kronig relations'^ 

e'M = l + -P r ^f^ds. (2) 

We assume that e"{uj) is narrowly peaked near a high (ultraviolet) freauency sq ^ w. 

At zero temperature e"{uj) may be a Dirac 5-function, or a sum of 5- functions, each 
corresponding to an allowed vibronic transition, or to its analogue in an amorphous solid. 
At finite temperature T the excited vibrational (phonon) substates of the ground electronic 
state are excited. In general, these have the same allowed transitions to excited electronic 
states as the vibrational ground state, and we will assume they have the same oscillator 
strengths. Assuming a classical distribution of ground state vibrational excitation (re- 
quiring, in a crystal, a temperature in excess of the Debye temperature) yields a simple 
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analytical model for e"{(jj): 
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(3) 



e"{s) = I ]£t exp[n(s - so)/kBT] s < sq 
[O s > So, 



where C is a normalization constant. This model assumes a single allowed electronic 
transition of (zero-temperature) frequency sq, but may be applied if there are many such 
transitions, each with a frequency ^ a;, if sq is taken as an oscillator strength- weighted 
average. It tacitly neglects any temperature dependence of the Franck-Condon factors. 
Using (3) in (2) yields 

2 s Ch 
e'{u) = 1 + -P / exp[h{s - so)/kBT] ds (4) 
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If we approximate kBT <^ Tisq then the integrand is narrowly peaked around s = sq. The 
lower limit of the integral may be extended to — oo and — a;^) expanded in a Taylor 

series about sq: 

+ o-:^(^-«o)'-.... (5) 
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This removes the pole in the integrand. Then the integral is elementary, yielding 

6'H-1 = ^^, (6) 
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which determines C. 

Differentiating (4) yields 

de'iuj) 2 r° s cn rn(s-so) \ , 

Again taking kBT -C ^sq and a; -C so, the first term in the expansion (5) integrates to 
zero but the second yields the result 

In terms of the refractive index n = (e//)^/^ = e^/^ 

dn{uj) — 1 kB 



dT 2n hso ' 
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(9) 



In order to evaluate this numerically we need to estimate sq. Fused silica has very strong 
ultraviolet absorption centered around a wavelength of 1200 A, implying Usq ~ 1.6 x 10~^^ 
erg. Then, with n = 1.445 in the visible and near-infrared, we have 

^ = 3.2 X 10-6 °K-^ (10) 

This theoretical value is about 3.4 times less than the measured^ value of 1.1 x 
10"^ °K-^. The most plausible explanation of the discrepancy is that the Franck-Condon 
factors are increasing functions of temperature, as increasing vibrational excitation of the 
electronic ground state increases the range of internuclear separations, providing a better 
match to upper states in which the equilibrium internuclear separation is greater. This 
effect cannot be calculated analytically. 

The theory predicts that dn/dT be nearly independent of frequency because n is nearly 
independent of frequency. This is observed^ to be true to within about 10% for wavelengths 
between 3000 A and 4//. dn/dT is measured to increase with decreasing wavelength for 
A < 3000 A, where the assumption a; <^ sq is no longer valid. The theory could readily be 
extended to these larger values of uj by means of a straightforward numerical integration. 
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